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ABSTRACT 

We report the discovery of an extremely curved jet in the radio-loud quasar PKS 2136+141. Multi- 
frequency Very Long Baseline Array (VLBA) images show a bending jet making a turn-around of 210 
degrees in the plane of the sky, which is, to our knowledge, the largest ever observed change in the 
position angle of an astrophysical jet. Images taken at six different frequencies, from 2.3 to 43 GHz, 
reveal a spiral-like trajectory, which is likely a sign of an intrinsic helical geometry. A space- VLBI 
image, taken with the HALCA satellite at 5 GHz and having comparable resolution to our ground- 
based 15 GHz data, confirms that the bend is a frequency-independent structure. VLBA monitoring 
data at 15 GHz, covering eight years of observations, show knots in the jet clearly deviating from 
ballistic motion, which suggests that the bending may be caused by a growing helical Kelvin-Helmholtz 
normal mode. The jet appearance suggests a helical wave at a frequency well below the "resonant" 
frequency of the jet, which indicates that the wave is driven by a periodic perturbation at the base 
of the jet. We fit the observed structure in the source with a helical twist, and we find that a simple 
isothermal model with a constant wave speed and wavelength gives a good fit. The measured apparent 
velocities indicate some degree of acceleration along the jet, which together with an observed change 
in the apparent half-opening angle of the jet allow us to estimate the changes in the angle between 
the local jet direction and our line of sight. We suggest that the jet in PKS 2136+141 is distorted by 
a helical Kelvin-Helmholtz normal mode externally driven into the jet (e.g. by precession), and that 
our line of sight falls within the opening angle of the helix cone. 
Subject headings: galaxies: jets — quasars: individual (PKS 2136+141) 



1. INTRODUCTION 

A significant fraction of extragalactic jets show some 
degree of bending - from slightly curved jets up to a 
complete turn- around of almost 180°. Recently, in their 
large s tudy of jet kinematics of r adio-loud active galactic 
nuclei, Kellermann et al. (2004) measured vector veloc- 
ities for 60 bright jet features (also called components). 
They found that approximately a third of these com- 
ponents show a significant non-radial motion, i.e. the 
direction of their velocity vector differs by at least 3ct 
from the mean structural position angle of the jet. If 
these observed velocities trace the underlying jet flow, 
their result indicates that bends in jet direction are very 
common. 

For core-dominated radio sources with high optical po- 
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larization, there is a well-known bimodal distribution of 
the angles between jets in parsec and kiloparsec scales, 
with a main peak of misalignm ent angles around 0° 
and a secondary peak aro und 90° l)Pearson & Readhead I 
119881 iLister et al. 1 120011) . How ever, a large-angle mis- 
alignment exceeding 120° is rare llWilkinson et al. 11198(1 
iTineav et al~l Il998t ILister et al I l2001fl . Up to today, 
the largest observed AP.A. is 177° in the gamma-ray 
blazar PKS 1510-089, which show s a jet bending almos t 
directly across our line of sight ijHoman el al. 1 1207)1 . 
Since core-dominated radio sources have jets oriented 
close to our line of sight, all intrinsic variations in the 
jet trajectories are exaggerated in projection - often to 
a large degree. This implies that rather small intrinsic 
bends can manifest themselves as large-angle misalign- 
ments between the jet axes observed on parsec and kilo- 
parsec scales, or as high as ~ 90° turns in the VLBI 
images. 

Observations of relativistic jets in parsec scales 
provide evidence that AG N jets can exhibit " wiggling" 
structures (e.g. 4C 73.18. IRoos et al. I Ip9% 3C345, 
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iZensus et all (11995ft: 3C27 3 lLobanov fc Zensusl 
<|2001ft : 3C120. lHardee et al. I l|2005ft ) reminiscent of 
helically twisted patterns. It has been proposed that 
Doppler boosting together with parsec-scale jets travel- 
ing in helical paths could explain the excess of sources 
showing 90° misalignment angle between pc and kpc 
scale jets wit hout invoking an uncomfor table 90° intrin- 
sic curvature l)Conwav fc Murphy 11993ft . As the number 
of sources showing apparently helical structures has 
grown, the helical jet models have become increasingly 
popular also as an explan ation for the (quasi)p e riodic 
flux variations in A GN l|Abraham fc Romero I 119991: 
lOstorero et al. 1 12004ft . However, the mechanism pro- 
ducing an apparently helical shape of the jet is unclear 

- as are the explanations also for more modest observed 
bends. 

The "corkscrew" structure of the jet in the well- 
known galactic source SS433 is successfully explained 
by ballistic motion of material ejected from a pro- 
cessing jet nozzle (|Stirling et al. I T2002ft . and a sim- 
ilar model has been suggested also for several ex- 
tragalactic jets show i ng " wiggling" ( e.g. 3C273, 
lAbraham fc Romero I <ll999ft: BLLac . iSthling et al. I 
(|2003ft : OJ287. lTatevama fc Kinghaml l)2004ft h In these 
models, the jet precession is either due to the Lense- 
Thirring effect in a case of misalignment between the 
ang ular momenta of accreti on disk and a Kerr black hole 
(see iCaproni et al. I l)2004ft and references therein), or 
due to a binary black hole system where a secondary 
black hole tidally induces the precession. 

Con trary to the above-mentioned cases, iLister et al. I 
(2003) report that although the powerful radio source 
4C +12.50 exhibits a jet ridge line highly reminiscent 
of that in SS433, it is most likely due to streaming in- 
stead of ballistic motion. Streaming helical motion arises 
naturally from spatial stability analysis of relativistic 
jets, since the jets are un stable against growing Kelvin- 
Helmholtz normal modes l)Hardee 11987ft . Provided there 
is a suitable perturbation mechanism present in the inner 
part of the jet, the distortion waves propagating down the 
jet can displace the whole jet (helical fundamental mode) 
or produce helically twisted patterns on the jet surface 
(fluting modes). If the jet carries a large scale electric 
current (so-called Poynting flux dominated jets), it is 
in addition unstable against magnetic kink instability, 
which could also produce o bserved "wiggling" structures 
l|Nakamura fc Meier ll200l . 

Not all jets with observed bends exhibit "wiggling" 
structure, and many of the observed changes in the jet 
direction can be explained without invoking helical mo- 
tions. Proposed explanations for curving jets include 
ram pressure due to winds in the intracluster medium, 
a density gradient in a transition to the intergalactic 
medium and deflections by massive clouds in the inter- 
stellar medium. Most likely, different mechanisms work 
in different sources. It would be valuable to be able 
to reliably identify the reason for bending in individ- 
ual sources, since the observed properties of the bend 

- correctly interpreted - can constrain several physical 
parameters of the jet and the external medium (see e.g. 
Hardee ( 200|) for the case of K-H instabilities). 

In this paper, we present Very Long Baseline Ar- 
ray (VLBA) images from a dedicated multi-frequency 



observation and from the VLBA 2 cm Survey 4 
l|Kellermann et al. 1 119981) showing that PKS 2136+141 
(OX161), a radio-loud quasar at moderately high red- 
shift of 2.427, has a parsec-scale jet, which appears to 
bend over 180° on the plane of the sky, being - to our 
knowledge - the largest ever observed change in the po- 
sition angle of an astrophysical jet (other sources show- 
ing very pronounced chang es in the jet dire ction include 
e.g. PKS 1510-089 lIHomari et al. I 12001 . 1803+784 
(jBritzen et al. 111999ft . and NRAO 150 (I. Agudo et al., 
in preparation)). In VLA images, PKS 2136+141 is 
a compact sourc e, showing no extend ed emission on 
arcsecond scales llMurphv et al. 111993ft. Both 5 GHz 
llFomalont et al. 112000ft and 15 GHz IjKellermann et al. I 
1998) VLBA observations reveal a core-dominated source 
with a short, sli ghtly bending jet. 

Originally, ITornikoski et al. I l|2001ft identified 
PKS 2136+141 as a candidate gigahertz-peaked spec- 
trum source (GPS) having a slightly inverted spectrum 
up to 8-10 GHz in the intermediate-to-quiescent state 
and a clearly inverted spectrum (a > +0.5) during 
outbursts. The high t urnover frequency reported in 
ITornikoski et al. I l|2001ft even puts PKS 2136+141 in 
the class of high f reque ncy peakers (HFP). Although 
iTorniainen et al.~l l)2005ft have recently classified the 
source as a fiat spectrum radio source having a convex 
spectrum only during outbursts, t he simultaneous con- 
tinuum spectra from RATAN-600 ijKovalev et al. 1119991 
S. Trushkin, private communication) do show a convex 
shape also in the intermediate-to-quiescent state, albeit 
with a little lower peak frequency. The source is variable 
at radio frequencies showing a factor of ~ 3 variations 
in cm-wavelength flux curves with characteristic time 
scale of ~ 5 — 6 years (see Figure ^| . The last strong 
outburst started around 1998 and peaked in late 2002 
and in early 2004 at 14.5 and 8 GHz, respectively. This 
indicates that both, our multi-frequency observations in 
2001 and 15 GHz VLBA monitoring during 1995-2004, 
caught the source during a major flare. 

The paper is organized as follows: the multi-frequency 
VLBA data demonstrating the 210° bend of the jet, to- 
gether with a space- VLBI observation from the HALCA 
satellite, are presented in § 2. In § 3, we present a 
kinematic analysis, derived from over eight years of the 
VLBA 2 cm Survey monitoring data at 15 GHz, indi- 
cating non-ballistic motion of the jet components. Also, 
changes in /3 a pp and in the apparent half-opening an- 
gle of the jet are investigated. In § 4, possible rea- 
sons for observed bending are discussed and a helical 
streaming model explaining the observed structure is pre- 
sented. Conclusions are summarized in § 5. Throughout 
the paper we use a contemporary cosmology with Hq = 
71kms~ 1 Mpc- 1 , Q M = 0.27 and n A = 0.73. For this 
cosmology and a redshift of 2.427, an angular distance of 
1 mas transforms to 8.2 pc and a proper motion of 0.1 
mas yr _1 to an apparent speed of 9.2 c. We choose the 
positive spectral index convention, S v oc v +a . 

2. MULTI-FREQUENCY OBSERVATIONS 

On May 2001 we made multi-frequency polarimet- 
ric VLBI observations of four HFP quasars, including 
PKS 2136+141, using the VLBA. Observations were split 

4 http:/ /www. cv.nrao.edu/2cmsurvey/ 
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Fig. 1. — Weekly averages of total flux density for PKS 2136+141 (OX 161) from UMRAO centimetre band monitoring. Observations 
at 14.5, 8.0 and 4.8 GHz are denoted by (green) crosses, (blue) circles and (red) triangles, respectively. Extrapolated ejection epochs of 
VLBI components arc marked with arrows. Component 16 is left out from the figure because it is seen only at the last epoch of the VLBA 
2 cm Survey monitoring data, and therefore, we cannot accurately determine its ejection epoch. The ejection of 16 must have occured at 
sometime between the ejection of 15 and April 2004, when it is first observed. See the on-line edition for a color version of the figure. 



TABLE 1 

Parameters of the Images for Figures [5] and [3] 



Frequency ®b,m&j ©i),min P.A. rms noise Peak Intensity Contour co a 

(GHz) (mas) (mas) (deg) (mjy beam -1 ) (mjy beam - 1 ) (mjy beam - ) 



2 


6.12 


3.02 


-11.7 


0.4 


1153 


1.2 


5 


3.37 


1.70 


-12.5 


0.2 


1811 


0.6 


5 b 


0.75 


0.59 


-26.6 


8.3 


708 


24.9 


8 


1.98 


1.01 


-7.8 


0.2 


1890 


0.6 


15 


0.90 


0.50 


-7.5 


0.5 


1429 


1.5 


22 


0.63 


0.36 


-7.9 


1.2 


1153 


3.6 


43 


0.43 


0.18 


-16.1 


1.5 


512 


4.5 



a Contour levels are represented by geometric series eo(l, 2 n ), where Co is the lowest contour 
level indicated in the table (3x rms noise). 
b HALCA image 



into a high frequency part (15, 22 and 43 GHz), which 
was observed on the 12th of May, and into a low fre- 
quency part (2.3, 5 and 8.4 GHz) observed on the 14th 
of May. Dual polarization was recorded at all frequen- 
cies. 

2.1. Reduction of the VLBA Data 

The data were correlated with the VLBA correlator in 
Socorro and were postprocessed in Tuorla Observatory 
using the N RAO's Astronomical Image Processing Sys- 
tem, AIPS, (|Bridle fc Greisen 1119941 iGreisen lll98cy and 
the Caltech difmap package l|Shenherd' Standard 



methods for VLBI data reduction and imaging were used. 
An a priori amplitude calibration was performed using 
measured system temperatures and gain curves. For the 
high frequency data (15-43 GHz), a correction for at- 
mospheric opacity was applied. After the removal of a 
parallactic angle phase, a single-band delay and phase 
off-sets were calculated manually by fringe fitting a short 
scan of data of a bright source. We did manual phase cal- 
ibration instead of using pulse-cal tones, because there 
were unexpected jumps in the phases of the pulse-cal 
tones during the observations. Global fringe fitting was 
performed, and the delay difference between right- and 
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Fig. 2. — Total intensity VLBA observations of PKS 2136+141 from May 2001. The compiled figure shows images of the source at all 
six observed frequencies (2.3, 5, 8.4, 15, 22, and 43 GHz). The letter "C" in the top left panel marks the location of the core. The size and 
orientation of the beam is shown in the lower left corner of each image. Peak intensities and contour levels are given in Table HI 



left-hand systems was removed (for the purpose of future 
polarization studies). Bandpass corrections were deter- 
mined and applied before averaging across the channels, 
after which the data were imported into difmap. 

In difmap the data were first phase self-calibrated us- 
ing a point source model and then averaged in time. 
We performed data editing in a station-based manner 
and ran several iterations of clean and phase self- 
calibration in Stokes I. After a reasonable fit to the clo- 
sure phases was obtained, we also performed amplitude 
self-calibration, first with a solution interval correspond- 
ing to the whole observation length. Solution interval 
was gradually shortened as the model improved by fur- 
ther cleaning. Final images were produced with the Perl 
library FITSPlot 5 . 

We have checked the absolute flux calibration by com- 
paring the extrapolated zero baseline flux density of our 
compact calibrator source 1749+096 at 5, 8.4, and 15 
GHz to the single-dish measurements made at University 
of Michigan Radio Astronomy Observatory (UMRAO), 
and at 22 and 43 GHz to the fluxes from Metsahovi Ra- 
dio Observatory's qu asar monitoring program at 22 and 
37 GHz, respectively ijTerasranta et al. 1 12004). The flux 
densities agree to 5% at 8.4, 15, 22 and 37/43 GHz, and 
to 8% at 5 GHz, which is better than the expected nom- 
inal accuracy of 10% for the a priori amplitude calibra- 
tion. Being unable to make a flux check for the 2.3 GHz 
data, we conservatively estimate it to have an absolute 

5 http : //personal . denison . edu/"homand/ 



flux calibration accurate to 10%. 

In order to estimate the parameters of the emission 
regions in the jet, we model-fitted to the self-calibrated 
(u,v) data in difmap. The data were fitted with a com- 
bination of elliptical and circular Gaussian components, 
and we sought to obtain the best possible fit to the visi- 
bilities and to the closure phases. Several starting points 
were tried in order to avoid a local minimum fit. We 
note that since the source structure is complex, the mod- 
els are not unique, but rather show one consistent pa- 
rameterization of the data. Based on the experiences 
in error estimation reported by several groups, we as- 
sume uncertainties in component flux ~ 5%, in position 
~ 1/5 of the beam size (or of the component size if it 
is larger than t he be am), and in size ~ 10% (see e.g. 
I.Torstad et a"l~l l|2005ft and iSavolainen et al. I l|2006|) for 
recent discussions on the model fitting errors). Although 
Uorstad et al. 1 1)2005|) use larger positional uncertainties 
for weak knots having flux densities below 50 mJy, we 
use ~ 1/5 of the beam size (or of the component size) 
also for these components. Bigger uncertainties would 
result in such a large ratio of the individual errors to 
the scatter of the component positions (about the best- 
fit polynomial describing the component motion) that it 
would be statistically unlikely (see § I3.1.2|l . 

A detailed description of the polarization data reduc- 
tion and imaging in Stokes Q and U together with po- 
larization images will appear in T. Savolainen et al. (in 
preparation). 
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Fig. 3.— 5 GHz space- VLBI image of PKS 2136+141 observed 
on 1998 May 28. Beam size, peak intensity and contour levels are 
given in Table HI 

2.2. Reduction of the HALCA data 

The 5 GHz space- VLBI observations of PKS 2136+141 
were carried out as a part of the VSOP survey program 
l)Hirabavashi i feOOO. R. Dodson et al. in preparation) on 
1998 May 28. In addition to the HALCA satellite, the ar- 
ray consisted of Arccibo, Haartebesthoek, and the Green 
Bank (140 ft) telescopes. The data were correlated us- 
ing the Penticton correlator and reduced using the same 
standard methods as for the ground based images. Due 
to the small aperture of the HALCA satellite, its weight 
was increased to 10 to persuade the fitting algorithm to 
take better account of the long space baselines. 

The dynamic range of the resulting image is rather 
small. This is due to a component that had a strong ef- 
fect only to a single baseline between Arecibo and Green 
Bank. Because the (u,v) coverage is very limited, all at- 
tempts to model this component diverged. By manually 
changing the model to follow the visibility amplitudes 
in this scan, we estimate that the position of this com- 
ponent is about 7 mas at a PA of 100°. Because this 
component could not be formally included in the model, 
the residuals are rather strong and thus the imaging noise 
is high. 

2.3. Source Structure from the Multi-Frequency Data 

Figure |21 displays images of PKS 2136+141 at all six 
observed frequencies. In the images at 15, 22 and 43 
GHz, uniformly weighted ( it, v)-grids are employed in or- 
der to achieve the best possible resolution, whereas nor- 
mal weighting is used in the low frequency maps to high- 
light diffuse, low surface brightness emission. The restor- 
ing beam sizes, peak intensities, off-source rms noise, and 
contour levels of the images are given in Tabled 

The multi-frequency images strikingly reveal a jet 
which gradually bends 210° with its structural PA. turn- 
ing clockwise from -27° at 43 GHz to +123° at 2.3 GHz. 
We identify the bright and the most compact model com- 
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Fig. 4. — Spectra of the core (filled circles) and the component 
15 (open squares) at 8-43 GHz on May 2001. 

ponent lying in the south-east end of the jet in the 15-43 
GHz images as the core, and mark it with a letter "C" . 
The identification is confirmed by a self-absorbed spec- 
trum of the component, as will be shown later. Between 
about 0.4-1.0 mas from the core, the jet turns over 90°, 
which is visible in 15-43 GHz images, and in the im- 
age taken at 8.4 GHz, the jet direction continues to turn 
clockwise ~ 50° at about 3.5 mas from the core. There 
is also evident bending in the 5 and 8.4 GHz images: a 
curve of ~ 70° takes place at about 6 mas south of the 
core. It is not totally clear, whether the trajectory of the 
jet is composed of a few distinct bends or whether it is 
a continuous helix. However, the gradual clockwise turn 
and the apparent spiral-like appearance of the jet in the 
multi-frequency images are highly reminiscent of a heli- 
cal trajectory. Also, it seems unlikely that the jet goes 
through at least three consecutive deflections with all of 
them having the same sense of rotation in the plane of 
the sky. 

The 5 GHz space- VLBI image of PKS 2136+141 from 
1998 May 28, i.e. three years before the multi-frequency 
VLB A observations, shows a rather compact core-jet 
structure with an extended emission to north-west (Fig- 
ure EJ. The jet takes then a sharp 90° bend to south- 
west within about 1 mas from the core. This 5 GHz 
image shows a very similar curved structure near the 
core that can be seen in our ground-based 15 GHz image 
with matching resolution. Hence, the observed large- 
angle bending between the images taken at different fre- 
quencies cannot be attributed to frequency dependent 
opacity effects. 

PKS 2136+141 is rather compact at all frequencies 
with a maximum jet extent of approximately 15 mas cor- 
responding to ~ 120 pc at the source distance. The two 
brightest model components are the core and a newly 
ejected component 15 (see §3 for component identifica- 
tions), which is located at ~ 0.25 mas from the core - 
near the beginning of the first strong bend in the jet. 
Figure 01 shows the 8-43 GHz spectra of these two com- 
ponents. At 2.3 and 5 GHz the angular resolution is 
too poor to separate the core from 15, and hence, corre- 
sponding flux values are omitted in the figure. The core 
has a synchrotron peak frequency of ~ 20 GHz, while 
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the newly ejected component 15 shows an optically thin 
synchrotron spectrum with a spectral index a = —0.6 
and a peak frequency below 8 GHz. The component 15 
is brighter than the core at every frequency except at 22 
GHz, where almost equal flux densities are measured. An 
optically thick spectrum of the core at frequencies below 
22 GHz indicates that 15 is brighter than the core also at 
2.3 and 5 GHz, although the spectra cannot be measured 
at those frequencies. The self-absorbed spectrum shown 
in Figure 01 confirms that we have correctly identified the 
core. 

3. VLBA 2 CM SURVEY AND MOJAVE DATA 

Seven observations of PKS 2136+141 have been 
performed as pa rt of the VLBA 2 cm Survey 
ijKellermann et al. I 119981 IZensus et all [2002) be- 
tween 1995 and 2004 of which the last one (epoch 
2004.27) 6 was part of t he follow -up program MOJAVE 
({Lister fc Homan"1 l2005). Details on the survey observ- 
ing strategy, the observing details and the data reduction 
can be found in the mentioned publications. 

clean VLBI images of PKS 2136+141 have been pro- 
duced by applying standard self-calibration procedures 
with DIFMAP for all seven epochs. The calibrated vis- 
ibility data were fitted in the (u,v) domain with two- 
dimensional elliptical Gaussian components (see Fig- 
ure EJ). Because of the complicated source structure, no 
adequate model representation of the source could be 
established with a smaller number of components than 
shown in Figure Moreover, the components are found 
to be located close-by along the curved inner part of the 
jet, in many cases separated from each other by consid- 
erably less than one beam size. Special care had to be 
taken to find a model of this complicated structure in a 
consistent way for all epochs. 

3.1. Source Kinematics 

A crucial question raised by the apparent helical shape 
of the jet in Figure [21 is, whether it represents streaming 
motion or whether the helix is due to a precessing jet 
nozzle ejecting material that moves along ballistic tra- 
jectories. We have plotted positions of all model-fit com- 
ponents from the multi-frequency data and from eight 
years of the VLBA 2 cm Survey monitoring data into 
FigureEl The open circles, which correspond to the mul- 
tiepoch data from the VLBA 2 cm Survey, form a dense 
and strongly curved region near the base of the jet, a ~ 2 
mas long continuous and slightly curved section 3.5 mas 
south-west of the core, and a few small isolated groups. 
This subtle finding alone, without invoking any compo- 
nent identification scenarios, demonstrates that ballistic- 
motion models are unlikely to yield a meaningful repre- 
sentation of the jet kinematics in PKS 2136+141. Bright- 
ening at certain points of the jet can be either due to an 
increased Doppler factor (if a section of the jet bends to- 
wards our line of sight), or due to an impulsive particle 
acceleration in a standing shock wave (e.g. forming in 
the bend) . There seems to be a zone of avoidance in the 
15 GHz data between the base of the jet and the large 
western group, further supporting this idea. 

6 After 2004.27, two further MOJAVE observations of 
PKS 2136+141 have been conducted. The monitoring will be con- 
tinued through 2006 and later. 



We have analyzed the jet kinematics using source mod- 
els derived from the VLBA 2 cm Survey monitoring data. 
In a complicated source like PKS 2136+141 it is often 
difficult to identify components across epochs with con- 
fidence. We have based our identifications on the most 
consistent trajectories, and on the flux density evolution 
of the components. We restrict our analysis to a full 
kinematical model for the inner 2 mas of the jet, since 
beyond that distance a fully self-consistent model could 
not be established due to the lower surface brightness 
and high complexity in the outer region. 

3.1.1. Component Trajectories and Flux Evolution 

Figure[7|shows the trajectories of the six compact com- 
ponents (11-16) which we have identified in the inner 2 
mas of the jet. The components travel towards south- 
west and their trajectories do not extrapolate back to 
the core implying that the components travel along a 
curved path, and the jet is non-ballistic. The compo- 
nents fade below the detection limit within ~ 1.6 mas 
from the core, and we cannot follow them further. 

We present the flux density evolution of the six inner- 
jet components in Figure |H1 The overall picture is that 
the fluxes decrease as the components move forward 
along their path. They reach about 5-15 mJy by the 
time they have traveled to a distance of ~ 1.5 mas from 
the core, after which they are not seen anymore. The 
knot 13 is the brightest component of the source at the 
first two epochs, and its flux density evolution matches 
well with the strong 1993 flare visible in the 14.5 GHz 
UMRAO flux density curve (see Figure^). The ejections 
of the components II and 12 (see Table 0J) take place 
during the rising phase of the 1993 outburst, and the 
component 13 is ejected just before the outburst reaches 
its maximum. 

The core has a steadily rising flux density until 2001.20, 
when the flux suddenly drops and two new components, 
14 and 15, appear. At the next two epochs, the core flux 
increases again and at the last epoch (2004.27) there is a 
drop accompanied by an appearance of a new component, 
16. The brightening of the core during our monitoring 
and ejections of components 14, 15, and 16 correspond 
to a strong total flux density flare peaking in late 2002 
at 14.5 GHz (again, see the total flux density curve in 
Figure P). 

The above-mentioned flux density evolution of the 
components is self-consistent and it is well in accordance 
with the general behavior of flat-spect rum radio quasars 
durin g strong total flux density flares ijSavolainen et al. 1 
2002), i.e. a compact VLBI core is mostly responsible for 
the rising part of the observed flares in single-dish flux 
curves, and a new component appears into the jet during 
or after the flare peaks accompanied by a simultaneous 
decrease in the core flux density. The fact that the over- 
all flux density evolution of the identified components 
in PKS 2136+141 seems to obey the common behavior 
identified for a number of other sources supports our 
kinematical model. An intriguing detail is the consec- 
utive ejection of, not one, but three new components in 
connection with a single, strong total flux density flare. 
This may indicate that there is some substructure in the 
flares, i.e. the outbursts in 1993 and 2002 could be com- 
posed of smaller flares. On the other hand, it could also 
mean that a single strong event in the total flux density 
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Fig. 5. — VLBA images of PKS 2136+141 at 15 GHz. The figure includes observations from the VLBA 2 cm Survey, from the MOJAVE 
Survey and a 15 GHz image from our multi-frequency data set observed in 2001.36. In all images, a uniformly weighted («,«)-grid is used. 
The Gaussian components fitted to the visibility data are shown as ellipses overlaid each image. The size and orientation of the beam is 
shown in the lower left corner of each image. Peak intensities and contour levels are given in Table |2*1 



TABLE 2 

Parameters of the Images for Figure [F] 



Epoch ®b,ma,j ©b.min P.A. rms noise Peak Intensity Contour co a 
(yr) (mas) (mas) (deg) (mjy beam -1 ) (mjy beam -1 ) (mjy beam -1 ) 



1995.57 


1.01 


0.53 


11.5 


0.9 


1254 


2.7 


1996.82 


0.92 


0.41 


-1.6 


0.8 


633 


2.4 


1998.18 


1.31 


0.61 


-3.4 


1.4 


929 


1.2 


1999.55 


1.45 


0.49 


-19.6 


0.7 


1429 


2.1 


2001.20 


1.14 


0.52 


-3.1 


0.4 


1153 


1.2 


2001.36 


0.90 


0.50 


-7.5 


0.5 


1429 


1.5 


2002.90 


1.15 


0.55 


-12.9 


0.4 


2120 


1.2 


2004.27 


1.09 


0.54 


-5.7 


0.3 


1745 


0.9 



a Contour levels are represented by geometric series co(l, ...,2 n ), where co is the lowest 
contour level indicated in the table (3x rms noise). 



curve is able to produce complicated structural changes 
in the jet, e.g. forward and reverse shocks could both be 
visible, or there could be trailing shocks forming in the 
wa ke of the main pertur bance as theoretically predicted 
by lAgudo et al. I <|2001f) and recently o bserved in sev- 
eral objects, most p rominently in 3C 111 ij.Torstad et al. I 
120051 iKadler II2005L Ka dler et al. in preparation) and in 
3C 120 l|G6mez et al. ll200l|) . 

3.1.2. Apparent Velocities and Acceleration 

Identification of the model-fit components across the 
epochs allows us to estimate the velocities of the compo- 
nents and to search for possible acceleration or decelera- 



tion. Since the components follow a highly curved path, 
we choose to measure the distance traveled along the jet 
by using an average ridge line, which is shown by the 
dashed line in Figure El The ridge line is determined 
by fitting a smooth function to the positions of model 
components within the inner 2 mas of the jet. 

Figure displays the motion of the components along 
the jet. To measure the component velocities and to esti- 
mate the ejection epochs, we have made standard linear 
least-squares-fits to the motions: 

l(U) = ao + a>i(U - t mid ), (1) 
where I is the traveled distance along the jet, ti is the 
epoch of observation (i — 1,...,N) and i m i<j = (t\ + 
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TABLE 3 

Best-Fit Polynomials for the Motion of Components 11-15 



Comp. 


k a 


«o 


«i 


(12 




u c 


11 


1 


1.50 ±0.10 


0.27 ±0.17 









12 


1 


1.32 ± 0.05 


0.17 ±0.02 




0.34 


5 




2 


1.29 ±0.08 


0.17 ±0.02 


0.01 ±0.01 


0.37 


i 


13 


1 


0.80 ± 0.04 


0.12 ± 0.02 




1.50 


6 




2 


0.68 ± 0.06 


0.13 ± 0.02 


0.02 ±0.01 


0.43 


5 


14 


1 


0.81 ± 0.05 


0.18 ± 0.04 




1.19 


2 


15 


1 


0.35 ± 0.06 


0.09 ± 0.04 




0.12 


2 



a Order of the fitted polynomial. 

b The x 2 value is relative; see text. 

c Number of the degrees of freedom in the fit. 
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Relative R.A. [mas] 

Fig. 6. — Model-fit component locations. Open circles are model 
components from the VLBA 2 cm Survey monitoring and filled 
symbols show the multi-frequency data: squares - 2.3 GHz, tri- 
angles - 5 GHz, diamonds — 8.4 GHz, circles - 15 GHz, stars - 22 
GHz, and asterisks - 43 GHz. The position of the core is marked 
with a cross. The inset is an enlarged plot of the innermost 2x1.5 
mas and it also displays an approximate ridge line of the jet. 

£at)/2. In order to search for a possible acceleration or a 
deceleration, we have also fitted the components 12 and 
13, which are observed over most epochs, by second-order 
polynomials: 

l(U)=ao + ai(U - t mid ) + a 2 (U - t mid ) 2 - (2) 

To assess the goodness of the fits, we have calculated a 
relative x 2 value for each fit. The x 2 values are relative, 
not absolute, since we do not strictly know the real la 
measurement errors of the component positions but have 
instead used 1/5 of the beam size projected onto the 
ridge line. Although the "1/5 of the beam size" estimate 
for the errors is rather well-justified in the literature, it 
is likely to be conservative. Thus the tabulated values 
of x 2 distribution do not provide a good statistical test 
of the goodness of the fits in our case. However, we 
assume that the beam sizes give a good estimate of the 
relative errors between the epochs, and hence, only one 
scaling factor (which is close to unity) for the positional 
errors (and for the \ 2 values) remains unknown. With 
this assumption, the relative \ 2 values can be used to 
compare the goodness of the fits against each other. The 
parameters of the fitted polynomials are given in Table 
together with relative reduced x 2 values and the number 
of the degrees of freedom, for each fit. 

We have gathered the average angular velocities along 
the jet, (fi) — a%, the average apparent speeds, (f3 app ), 
and the epochs of zero separation, To, for each compo- 
nent in Table 0] For the second-order fits, also the an- 
gular acceleration along the jet, fi = 2a2, is reported. 
The proper motions range from 0.09 to 0.27 mas yr" 1 , 
corresponding to the apparent superluminal velocities of 
8.7 -25.1c. 

While the component 12 is well-fitted by a straight line, 
the linear proper motion model does not seem to ade- 
quately represent the motion of the component 13 as can 
be seen from FigureEl On the other hand, a second-order 
polynomial gives an acceptable fit to 13, and its relative 
\ 2 value divided by the number of the degrees of freedom 



is by a factor of 3.5 smaller than the value for the first- 
order polynomial. However, the small number of data 
points and the uncertainty about the real (lc) errors in 
the component positions makes it difficult to give a sta- 
tistical significance of the deviation from the constant 
speed model. We have used two approaches to assess the 
significance. The most robust and the most straightfor- 
ward test is to use a statistic F Vl Va = (Xi/^i)/(xi/ 1/ 2), 
which has an ^-distribution with (v\, v-i) degrees of free- 
dom, and which is independent of the unknown scaling 
factor in the positional errors. This statistics can be 
used to test whether the squared residuals of the linear 
model are significantly larger than those of the acceler- 
ating model. As already mentioned, for the first and the 
second-order polynomials in the case of 13, = 3.5, 
which implies a difference in the residuals only at the 
confidence level of a = 0.10; i.e. according to the F-test, 
the statistical significance of the difference in the resid- 
uals is only marginal. Another approach to this prob- 
lem is to try to estimate the real la uncertainties of the 
component positions by applying a method described by 
iHoman et al. 1 1)2001[) . They estimated the uncertainties 
of the fitted parameters of their proper motion models by 
using the variance about the best-fit model, and as a by- 
product they obtained an upper bound estimate of the 
component position uncertainty Since the relative % 2 
values divided by the degrees of freedom are significantly 
less than unity for components 12 and 15 (the situation 
of component 14 is discussed later), we suspect that the 
real la positional uncertainties are in fact smaller than 
1/5 of the beam size. If the beam sizes give a good es- 
timate of the relative errors between the epochs, we can 
try to estimate the uniform scaling factor for the po- 
sitional uncertainties by requiring that x 2 — v f° r the 
linear proper motion models of 12, 14, and 15. The sum 
T,i(x 2 /k 2 — Vi) 2 , (i = 12, 14, 15) is minimized for the posi- 
tional error scaling factor k = 0.77; i.e. the la positional 
uncertainties are 0.77x 1/5 of the beam size ss 0.08 — 0.15 
mas. If these (la) uncertainties are assumed, the linear 
fit to the motion of component 13 has x 2 = 15-2 with 6 
degrees of freedom, implying that the speed of 13 is not 
constant at the significance level of a — 0.025. 

Component 14 seems to show motion that indicates 
acceleration similar to 13, but since there are only four 
data points for 14, we have not fitted it with a second- 
order polynomial, which would leave only one degree of 
freedom for the fit. By looking at Figure one can see 
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Fig. 7. — Trajectories of components 11—16. Positions are measured relative to the core, which is placed into the origin in each panel 
(filled circle). The errorbars correspond to 1/5 of the beam size. 
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Epoch [year] 

Fig. 8. — Flux densities of the identified components from the 
VLBA 2 cm Survey monitoring. 



that 14 should have been present in the jet already at 
the epoch 1999.55, but it cannot be unambiguously in- 
serted into the source model; i.e. adding one more com- 
ponent in the model does not improve the fit, but rather 
makes the final fit highly dependent on the chosen ini- 
tial model parameters. One possible explanation is that 
the proper motion for 14 is faster than indicated by the 
best-fit line in Figure |3 and it is not actually present 
in the jet in 1999.55. A second possibility is that 14 is 
a trailing shock released in the wake of the primary su- 
pcrlumina l component 13 instead of being ejecte d from 
the core ijAeudo et al. II200H iGomez et al. II2001|) . This 
alternative is n ot a v ery plausible one, since according to 
lAgudo et al. I 1)200 1|) trailing shocks should have a sig- 
nificantly slower speed than the main perturbance. A 




1995 1996 1997 1998 1999 2000 2001 2002 2003 2004 2005 
Epoch ]year] 

Fig. 9. — Component motion along the jet ridge line. The 
symbols for different component are: open circles - II, filled squares 
— 12, open diamonds - 13, crosses - 14, filled triangles - 15, and an 
open square - 16. The solid lines represent linear fits to the data, 
and the dashed line represents a second-order fit. 

third possible explanation for the non-detection of 14 in 
1999.55 is that either 14 and the core or 14 and 13 are so 
close to one another in 1999.55 that they do not appear 
as separate components in our data. 

As discussed above, we cannot firmly establish that 
component 13 accelerates along the jet ridge line, al- 
though it seems probable. However, also the other com- 
ponents (apart from 14, which itself might show accelera- 
tion) appear to be systematically faster, the farther away 
from the core they lie. We have averaged the velocities 
of the individual components over traveled distance bins 
of 1.0 and 0.5 mas. In averaging, we have taken into ac- 
count the standard errors given in Table 0] and the fact 
that for a given component in a given bin, the weight 
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TABLE 4 

Proper Motions of Model-Fit Components 



Comp. 


fc a 


A 


0*) 




To 






(mas yr -2 ) 


(mas yr 


(c) 


(yr) 


11 


1 




0.27 ±0.17 


25.1 ± 15.2 


1990.7 ±3.0 


12 


1 




0.17 ±0.02 


15.8 ± 1.8 


1991.5 ±0.6 




2 


0.01 ±0.02 


0.17 ±0.02 


16.0 ± 1.8 




13 


1 




0.12 ±0.02 


11.2 ± 1.9 


1993.3 ±0.6 




2 


0.03 ±0.01 


0.13 ±0.02 


11.7 ±1.4 




14 


1 




0.18 ± 0.04 


16.2 ±3.9 


1998.2 ±0.8 


15 


1 




0.09 ± 0.04 


8.7 ±4.1 


1999.0 ±1.2 




0.4 0.6 0.8 1 1.2 1.4 
Projected distance along the jet [mas] 

Fig. 10. — Apparent velocity of the components as a function of 
projected distance along the jet. The figure shows two different bin 
widths: open squares — a bin width of 1.0 mas, and filled circles - 
a bin width of 0.5 mas. 



has to be decreased, if the first or the last observation of 
the component falls into that bin. This is because natu- 
rally we do not have any knowledge about the component 
speed before or after our observations. For component 
13 we have used the quadratic fit instead of the linear, 
since it much better describes the data. The binned ve- 
locities are presented in Figure ITU1 which suggests that 
there is acceleration along the jet ridge line. When the 
data is divided in four bins, the first bin (0-0.5 mas) has 
an average component velocity of ss 7±3 c, while the rest 
of the bins (0.5-2.0 mas) have significantly higher veloc- 
ities, the maximum being « 17 ± 3 c for the bin covering 
1.0-1.5 mas. For the division in two bins, the change 
is from sw 11 ± 1 c in the first bin to « 17 ± 3 c in the 
second. The acceleration takes place in the first strong 
bend of the jet, which seems very natural since the ap- 
parent velocity depends on the angle between the local 
direction of the jet flow and our line of sight. The accel- 
eration in Figure HTH can alternatively be explained with 
the component speed changing as a function of ejection 
epoch - a new component (15) being slower than the old 
ones (11-12). This, however, does not explain the likely 
acceleration of 13. Thus we regard the change in /3 app 
along the jet ridge line as probable although more data 



s 15 



Note. — See text for the definitions of ft, {fi), (/3 a pp), and To. 
a Order of the fitted polynomial. 
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Projected distance along the jet fmas] 

Fig. 11. — Apparent half-opening angle >J7app of the components 
as a function of projected distance along the jet. The symbols for 
components 11-16 are the same as used in Figs. 6 and 7. Stars 
correspond to components of the 22 GHz data and asterisks are 
for the 43 GHz components. The dashed line shows the average 
^Papp in two regions: 1) the innermost 0.5 mas, and 2) between 
0.8-2.0 mas. 

is needed for conclus i ve evid ence. 

iKellermann et al~l ((2004) parameterized the 1995 — 
2001 VLBA 2 cm Survey data of PKS 2136+141 by fit- 
ting Gaussian components to the brightest features of 
each epoch in the image plane and derived a speed of 
Pa,pp — 1-8 ± 1.4, which is much lower than our esti- 
mates. Basically, their approach traces the component 
13 over the time when its proper motion was slow (see 
Figure 13, and due to the curved path of 13, the speed 
measured from the changes in its radial distance from 
the core is also lower than its true angular velocity. A 
linear least-squares-fit to the radial distances of 13 in our 
data over the epochs 1995 — 2001 yields /3 an n = 2.6 , which 
agrees with the speed given in IKellermann et a l. (2004) 
within the errors. 

3.2. Apparent Half-Opening Angles and Jet Inclination 

The angle between the local jet direction and our line 
of sight affects also the apparent opening angle of the jet. 
Assuming that the jet has a conical structure and a filling 
factor of unity, we have estimated the local value of the 
apparent half-opening angle for component i as \]/^ pp = 
tan _1 (tij/Zj), where U is the distance of the component 
i from the core measured along the jet ridge line, and 
di = ai/ 2 + Ar i5 with <Zj being the size of the elliptical 
component projected on the line perpendicular to the 
jet and Ar^ being the normal distance of the component 
from the jet ridge line. 

Figure 1111 shows the apparent half-opening angle for 
components 11-16 from 15 GHz data and for unnamed 
components from 22 and 43 GHz data as a function of 
traveled distance along the ridge line. The components 
with zero axial ratio have been excluded from the figure. 
Naturally, 'J'app is very uncertain for a single component, 
but Figure ^2 clearly shows that there is a change in the 
average ^app between the first ~ 0.5 mas and 0.8 — 2.0 
mas, and the change seems to correspond to the first 
large bend in the jet. Within the first 0.5 mas, the av- 
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erage apparent half-opening angle is ~ 25° and further 
down the jet it decreases to ~ 10°. A natural explanation 
for this effect is that the viewing angle 9 increases in the 
first apparent bend, since ^app ~ *int / sin 9 if 9 and the 
intrinsic half-opening angle ^int are both small. There is 
a caveat, though. Namely, the large value of ^app near 
the core could be due to a resolution effect: if the po- 
sitional uncertainties for the components near the core 
were ~ 1/5 of the beam size, Ar of these components 
could be « 0.1 — 0.2 mas at 15 GHz merely due to the 
scatter, and this would result in (\l/a PP ) = 25°. However, 
there are two reasons, why we think the observed change 
in \& apP is a true effect and due to a changing 9. First, 
the model-fit components from 22 and 43 GHz data show 
similar values for \& app near the core as the components 
from 15 GHz data, although they have a factor of 1.5-2 
smaller positional errors. Secondly, the components near 
the core are bright, having fluxes comparable to the core, 
which makes their positional uncertainties smaller than 
the adopted ~ 1/5 of the beam size (see also the previous 
section for a discussi on about the position al uncertain- 
ties). According to Uorstad et al. I l)2005|) bright (flux 
of the knot > 100 rms noise level) and compact (size 
< 0.1 mas) features in the jet have positional uncertain- 
ties ~ 0.01 mas at 7 mm. Thus we regard the observed 
change in ^app as a genuine effect. 

The f3 app of the jet seems to increase after ~ 0.5 mas 
from the core (see Figure HD)l while the ^app decreases. 
This can be understood, if the angle between local jet 
direction and our line of sight within the first ~ 0.5 mas 
is smaller than the angle 9sl = sin (1/F), which max- 
imizes /3 a pp- After 0.5 mas 9 increases towards the max- 
imal superluminal angle increasing /3 a pp and decreasing 
V^app- Assuming that the largest velocity in Figure E3 
Papp ~ 17 c, is close to the maximum apparent velocity 
(we do not consider the velocity of II, 25.1 c, very reliable 
since it is based on only two data points, and hence, we 
do not use it in our estimation of the maximum /3 app ), 
we can estimate that the jet Lorentz factor T ~ 20. Now, 
if we assume a constant Lorentz factor T = 20, we can 
fit for the values of the viewing angle and the intrin- 
sic jet half-opening angle giving the observed superlumi- 
nal speeds and half-opening angles before and after the 
bend. For (/3 apP ) ~ 7 c, (* ap p) ~ 25° within the first 
0.5 mas from the core and (/3 app ) ~ 17 c, (^app) ~ 10° 
after the bend, we obtain the following best-fit values: 
(9(0 - 0.5 mas) = 0.6°, 9(1.0 - 2.0 mas) = 1.5°, and 
^int = 0.26°. According to this result, the jet bends 
away from our line sight by « 0.9° after the first 0.5 mas 
from the core. 

Naturally, it is possible that the components exhibit 
also real acceleration along the jet in addition to the 
apparent acceleration due to the projection effect. Un- 
fortunately, the observations do not allow us to decide 
between these cases. Applying Occam's razor, we con- 
sider a constant jet speed in the following discussion. 

4. DISCUSSION 

4.1. Possible Reasons for Bending 

In this section we discuss the mechanisms capable of 
producing the curved structure of PKS 2136+141. In 
principle, there are several possible scenarios for the ob- 
served bending in relativistic jets, but we can rule out 



some of them in this particular case on the grounds of 
our analyzed data. 

Projection effects play an important role in this case 
making the intrinsic bending angle much less than the 
observed one. If the angle between the jet and our line 
of sight is as small as the analysis in the previous section 
indicates, it is possible that, for example, within 0.5-1.0 
mas from the core where the apparent bend is ~ 90°, 
the intrinsic bending angle is only ~ 1 — 2 degrees. This 
implies that formation of internal shocks in the bend, and 
their possibly destructive effect to the collimation of the 
j et, is not of great concern here. iMendoza fc Longair I 
(2002) calculate an upper limit to the bending angle of 
a jet in order not to create a shock wave at the end 
of the curvature. Their result for relativistic jets is ~ 
50°, which leaves our estimated intrinsic bending angle 
for PKS 2136+141 well within the limits. Therefore, we 
do not consider internal shocks formed by bending to 
restrict the possible explanations for the curvature in this 
source. 

One of the scenarios we can rule out, is a precessing jet 
where the components ejected at different times to differ- 
ent directions move ballistically, and form an apparently 
curved locus. Such models have been used to describe 
oscillating 'nozzles' observed in some BL Lac sou rces 
{Stirling et al. 1120031 iTatevama fc Kingham 1 120041) . In 
the case of PKS 2136+141, we can reject the precess- 
ing ballistic jet hypothesis, since the individual compo- 
nents do not follow ballistic trajectories, but rather ex- 
hibit streaming motion along a curved path (see fl3.1.1|) . 
However, precession of the jet inlet may still be behind 
the observed structure if it serves as an initial pertur- 
bation driving a helical Kelvin-Hclmholtz normal mode 
( seej^SL 

iHoman et al. I ((2002) explained the large misalignment 
between the pc and kpc scale jets of PKS 1510-089 with 
a scenario where the jet is bent after it departs the 
host galaxy, either by a density gradient in the tran- 
sition region or by ram pressure due to the winds in 
the intracluster medium. In PKS 2136+141 the bend- 
ing starts within 0.5 mas from the core, and because 
9 = 0.6+oJ degrees for the inner part of the jet (the 
error range here refers to an uncertainty introduced by 
errors in (/3 app ); see § 13.1.21 and I3.2JI . the corresponding 
deprojected distance is smaller than 0.8 kpc. The whole 
bending visible in Figure takes place within about 15 
mas from the core. Given the viewing angles estimated 
in § 13.21 it is fair to say that the deprojected length 
of the jet is - probably significantly - less than ~ 12 
kpc, and the whole observed bending takes place within 
that distance from the core. There are no observations 
of the host galaxy of PKS 2136+141, but the depro- 
jected lengths estimated above can be compared with 
typical scalelengths of the elliptical hosts of radio-loud 
quasars (RLQs), which are reported by seyeral g roups. 
In their near-infrared study, iTavlor et al. I l)1996fl found 
low-redshift (z ~ 0.2) RLQ hosts to have very large half- 
light scalelengths R1/2', from 14.0 to 106.8 kpc with an 
average of ~ 30 kpc. Other stu dies have yielded smaller 
values: e.g. iFlovd et al. I <|2004|) used Hubble Space Tele- 
scope WFPC2 data to study hosts of 17 quasars at 
z ~ 0.4 finding (R 1/2 ) = 10.2 ±1.8 kpc for RLQs, and 
iKotilainen et al. I l)1998f ) reported (R1/2) = 13±7 kpc for 
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12 flat spectrum radio quasars up to z = 1.0 in their near- 
infrared study. In the framework of hierarchical models 
of galaxy formation, the hosts of quasars at high red- 
shift are also expected to be more compact th an their 
low-redshift counterparts. iFalomo et al. I l)2005fl have re- 
cently managed to resolve a quasar host at z = 2.555 
and they report an effective radius of 7.5 ±3 kpc. As the 
deprojected distance < 0.8 kpc indicates, the bending 
of the jet in PKS 2136+141 starts well within the host 
galaxy, and the typical scalelengths listed above suggest 
that most of the curved jet is located within the host, al- 
though it is possible that part of it lies in the outskirts of 
the host galaxy being susceptible to the density gradient 
in the transition region. However, the bending clearly 
starts inside the galaxy, and thus requires some other 
cause. 

A collision inside the host galaxy, between the jet and 
a cloud of interstellar matter, can change the direction of 
the flow. For example, iHoman et al. I l|2003() have found 
component C4 in 3C 279 to change its trajectory by 26° 
in the plane of the sky, and they suggest this to be due to 
a collimation event resulting from an interaction of the 
component with the boundary between the outflow and 
the interstellar medium. However, the mere fact that 
the observed AP.A. is larger than 180° in our case con- 
strains the collision scenario since it excludes the situa- 
tions where the jet is bent by a single deflection, e.g. from 
a single massive cloud. In principle, the observed struc- 
ture may be due to several consecutive deflections from 
a number of clouds, but, as already mentioned in H2.3I 
this would require there to be at least three successive 
deflections having the same sense of rotation in the plane 
of the sky. This is unlikely, although not impossible. Fu- 
ture observations aimed at detecting emission from the 
jet beyond 20 mas will be very interesting, since they 
could tell whether the bending continues in the larger 
scale or not. If curvature with the same sense of rotation 
as within 15 mas is observed also in the larger scale, the 
scenario with multiple deflections from ISM clouds seems 
highly unlikely. 

4.2. Helical Streaming Model 

Rclativistic jets are known to be Kelvin-Helmholtz un- 
stable, and they can naturally develop helical distortions 
if an initial seed perturbation is present at the jet ori- 
gin. The helical K-H fundamental mode is capable of dis- 
placing the entire jet and consequently producing large 
scale helical st ructures where the plasma streams along 
a curved path ijHardee Ill987|) . The initial perturbations 
can have a random spectrum, and be due to for exam- 
ple jet-cloud interaction, or they can originate from a 
periodic variation in the flow direction of the central 
engine (precession or orbital motion). These perturba- 
tions can trigger pinch, helical or higher order normal 
modes propagating down the jet, and the appearance of 
this structure depends on the original wave frequency 
and amplitude, as well as subsequent growth or damp- 
ing of the modes. Since the individual components in 
PKS 2136+141 do not follow ballistic trajectories, but 
rather seem to stream along a helical path, we consider 
a helical K-H fundamental mode to be a possible expla- 
nation for the curious appearance of this source. For 
simplicity, in the following discussion we limit ourselves 
to a purely hydrodynamical case and do not consider the 



effect of magnetic fields for the growth of K-H modes, 
nor discuss the current-driven instabilities, which can 
produce helical patterns in Poynting-flux dominated jets 
ijNakamura fc Meier ll2004]) . 

The appearance of the jet in PKS 2136+141 already 
suggests some properties of the wave. The fact that com- 
ponents are observed to follow a nearly stationary helix 
implies that the wave frequency of the helical twist, u, 
is much below the "resonant" (or maximally unstable) 
frequency w*, which corresp onds to the freq uency of the 
fastest growing helical wave (|Hardee 1119871) . Such a low 
frequency, long wavelength helical twist suggests that the 
wave is driven by a periodic perturbation at the jet base. 
If the central source induced white-noise-like perturba- 
tions, we would expect to see structures corresponding 
to th e fastest growing frequency, i.e. lo* (H ardee et al. I 
1994). 

If isothermal jet expansion without gradients in the 
jet speed or in the ratio of the jet and external medium 
sound speeds is assumed, the wave speed, f3 w (in units of 
c), in the low frequency limit remains constant as the jet 
expands (Hardee 2003). Since the intrinsic wavelength, 
A, for a given u> varies proportional to the the wave speed, 
it can be also assumed to remain constant along the jet 
(as long as ix> <C uj*). Applying this assumption, we have 
fitted a simple helical twist to the observed positions of 
the VLBI components. The helical twist is specified in 
cylindrical coordinates with z along the axis of helix, by 
an amplitude A in the radial direction, and by a phase 
angle cf> given by 

r z d z 

4> = 2-Kh TT^+01, (3) 
Jzi A \ z ) 

where h is the handedness of the helix (—1 for right- 
handed and +1 for left-handed), and </>i = 4>(z\). With 
constant A, the integral in equation Q becomes trivial 
and (j> oc z. The amplitude growth is assumed to be 
conical: A — ztantpc, where ^ c is the opening angle of 
the helix cone. To describe the orientation of the helix in 
the sky we use two angles: i is the angle between the axis 
of the helix and our line of sight and x is the position 
angle of the axis of the helix projected on the plane of the 
sky. There are altogether five parameters in the model 
since the handedness can be fixed to h = +1 by simply 
looking at the images. 

A computer program was written to fit the positions of 
the VLBI components with different models expected to 
explain the shape of the jet and the kinematics of individ- 
ual components. In the program, the cumulative sum of 
perpendicular projected image-plane distances between 
a three-dimensional model and the component centroids 
from the VLBI ob servations was minim ized using a dif- 
ferential evolution l)Storn &: Price Il997j) algorithm. This 
evolutionary algorithm was chosen because of its good 
performance with non-linear and multimodal problems. 
All the observed components (and all frequencies) were 
used in a single model- fit, because the components seem 
to follow a common path, i.e. the helical structure ap- 
pears to persist for longer than the component propaga- 
tion times. In addition to component positions on the 
plane of the sky, the model was constrained by the view- 
ing angles determined in M3.2I The angle 8 between the 
normal to the jet's cross section and our line-of-sight was 
evaluated along the model, and it was required to be com- 
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Fig. 12. — Best-fit trajectory and component centroid positions at different frequencies for the model describing a low-frequency helical 
fundamental mode of Kelvin-Helmholtz instability. 



TABLE 5 

Best-fit Parameters of the Helical Streaming Model 



Parameter 




Symbol 


Value 


Initial helical wavelength [mas] 




Ai 


776 


Half-opening angle of the helix 


cone [deg 


lj>c 


1.0 


Initial phase angle [deg] 






-69 


Viewing angle to cone axis [deg 


] 


i 


0.3 


Sky position angle of cone axis 


[deg] 


X 


-158 


Helix handedness 




h 


+1 



patible with the viewing angles determined in M3.2I The 
following limits for 9 were used: at 0.3 mas from the core 
9 is between 0.4 — 0.8°, and at 1.5 mas from the core 9 is 
between 1.3 — 1.7°. If the model exceeded these limits, 
the cost function in the fitting algorithm was multiplied 
with a second-order function normalized to the width of 
the allowed range. 

The best-fit trajectory of the model is presented in Fig- 
ure El an d Figure ITTA shows the angle between the local 
jet direction and our line of sight, the Doppler factor, and 
the apparent velocity as a function of distance along the 
jet, which have been calculated by assuming T — 20. As 



is clear from these figures, the simple helical twist with 
constant wavelength along the jet gives a good fit to the 
data. The best-fit parameters of the model are listed in 
Table [5] In the best-fit model, we are looking straight 
into the cone of the helix (i = 0.3°), which has a small 
half-opening angle of 1.0°, meaning that the orientation 
of the helix is a very lucky coincidence. The fitted helical 
wavelength of the perturbation is 776 mas, corresponding 
to 6.4 kpc. This still needs to be corrected for a combined 
effect of the geometry and the possibly relativistic wave 
speed. The true intrinsic helical wavelength is given by 
Aint = A mo doi(l — f3 w cosi). Some constraint for (3 W can 
be derived from the fact that we do not see any system- 
atic change in the helical trajectory during 8.7 years, i.e. 
the change of the trajectory is less than the components' 
positional uncertainty in our observations. This gives an 
upper limit on the apparent wave speed, /32P P < 1.1c. 
As we look into the helix cone, the appropriate viewing 
angle for the wave motion is between ip c — i and i\) c + i, 
which gives (3 W < 0.989 and A int > 0.01A mo dci- These 
limits are purely due to the uncertainty in the positions 
of the VLBI components (~ 1/5 of the beam size) and 
it is likely that the true wave speed is much slower and 
the intrinsic wavelength is closer to A mo doi- Further ob- 
servations, e.g. within MOJAVE program, should give a 
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Fig. 13. — Viewing angle (top), Doppler factor (second from 
the top), apparent component speed (third from the top) and de- 
projected distance (bottom) for a jet with T = 20 as a function of 
pro ject ed distance from the core along a helical model (see Fig- 
ure 1121 . The two small vertical bars in the top panel show the 
viewing angle constraints obtained from the source kinematics and 
used in the modeling. 



tighter constraint for (3 W . In the low frequency limit, the 
wave speed is 

T 2 n 



i + rv 



where f3j = yl — T 2 is the flow velocity, and rj = 
(a x /ctj) i , with Oj and a x being the sound speeds in the je t 
and in the external medium, respectively ijHardee 1 2003). 
Assuming V = 20 and applying the upper limit of (3 W , we 
get a limit on the ratio of sound speeds: a x /aj < 0.5. 

We would like to have an estimate of the wave fre- 
quency of the observed helical twist since it could pos- 
sibly tell us about the origin of the periodic perturba- 
tion. Unfortunately, we do not know the wave speed, 
and hence cannot calculate the frequency. The lim- 
iting case with (3 W = 0.989 considered above yields 
u ~ 1.4 x 10~ 10 Hz corresponding to a period of ~ 200 yr, 
which is likely much below the actual period. However, 
some example values can be calculated for different com- 
binations of a x and Oj. Let us first assume an ultra- 

relativistic jet with ctj = c/v3 and F — 20. Now, using 
equation (0} we can calculate (3 W an d consequently oj for 
differe nt values of a x . For example, IConwav fc WrobeTl 
l|1995|) use a x » 400 km/s in their study of Mrk501 
on basis of X-ray observations and theoretical model- 



ing of giant elliptical galaxies, but this value refers to 
an average sound speed in the central regions of an el- 
liptical galaxy. Around the relativistic jet, there may 
be a hot wind or cocoon, where the sound speed is 
much higher, being a significant f ractio n of the light 
speed. For instance. lHardee et a l. (2005) estimate that 



the sound speed immediately outside the jet in the ra- 
dio galaxy 3C 120 is a x > 0.1c. For these two values, 
a x = 0.001c and a x = 0.1c, we get to ~ 2 x 10~ 15 Hz and 
lo r~j 2 x 10~ 13 Hz, respectively. The corresponding driv- 
ing periods are Pd ~ 10 7 yr and Pd ~ 10 5 yr. If ctj is less 

than c/-\/3 of the ultra-relativistic case, the frequencies 
will be higher and corresponding periods shorter. 

The helical streaming model presented above describes 
only one normal mode and it does not explain why the 
components in the 15 GHz monitoring data cluster in 
groups, i.e. why there seem to be certain parts in the jet 
where the flow becomes visible (see § 3.1). However, this 
might be explained if there are other, short wavelength, 
instability modes present in the jet in addition to the 
externally driven mode. Our current data set does not 
allow to test this hypothesis, but some hints of another 
instability mode may be present in Figure IT21 where the 
components look like they are "oscillating" about the 
best-fit trajectory. 

5. CONCLUSIONS 

We have presented multi-frequency VLBI data reveal- 
ing a strongly curved jet in the gigahertz-peaked spec- 
trum quasar PKS 2136+141. The observations show a 
210° change of the jet position angle, which is, to our 
knowledge, the largest ever observed AP.A. in an astro- 
physical jet. The jet appearance is highly reminiscent of 
a helix with the axis of the helix cone oriented towards 
our line of sight . 

Eight years of monitoring with the VLBA at 15 GHz 
show several components moving down the jet with 
clearly non-ballistic trajectories, which excludes the pro- 
cessing ballistic jet model from the list of possible expla- 
nations for the helical structure in PKS 2136+141. In- 
stead, the individual components are streaming along 
a curved trajectory. The estimated ejection epochs of 
the components are coincident with two major total flux 
density outbursts in 1990's, with three components be- 
ing associated with both outbursts. This may suggest 
that a single flare-event is associated with complicated 
structural changes in the jet, possibly involving multiple 
shocks. 

Most of the components have apparent velocities in the 
range of 8 — 17 c. One component even shows 25 c, but 
this high speed is based on only two data points, and 
therefore it is unreliable. We find evidence suggesting 
that /3 a pp increases after the first 0.5 mas from the core, 
a distance corresponding to a strong bend. Since the 
apparent jet opening angle also changes at this point, we 
suggest that the angle between the local jet direction and 
our line of sight increases from ss 0.6° mas within first 
0.5 mas to w 1.5° at distances between 1 and 2 mas. 

We fit the observed jet trajectory with a model de- 
scribing a jet displaced by a helical fundamental mode 
of Kelvin-Helmholtz instability. Our observations sug- 
gest that the wave is at the low frequency limit (relative 
to the "resonant" frequency of the jet) and has a nearly 
constant wave speed and wavelength along the jet. This 
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favors a periodic perturbation driven into the jet. The 
source of the perturbation could be e.g. jet precession 
or orbital motion of a supermassive binary black hole. 
Our present data does not allow us to reliably calculate 
the driving period of the K-H wave, and hence we can- 
not further discuss the perturbation's origin. Follow-up 
observations of the jet at angular scales larger than 15 
mas, as well as further monitoring of the jet kinematics, 
will probably shed light on this question. In the best- 
fit model, the helix lies on the surface of the cone with 
a half-opening angle of 1.0°, and the angle between our 
line of sight and the axis of the cone is only 0.3°, i.e. we 
are looking right into the helix cone. 
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